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Abstract 

The multi-electron wave function of an interacting electron system depends on the 
size of the system, i.e. the number of electrons. Here the question investigated is 
how the wave function changes for a symmetric Friedel- Anderson impurity when the 
volume is doubled. It turns out that for sufficiently large volume (when the level 
spacing is smaller than the resonance width) the change in the wave function can be 
expressed in terms of a universal single-electron state \q) centered at the Fermi level. 
This electron state is independent of the number of electrons and independent of the 
parameters of the Friedel- Anderson impurity. It is even the same universal state for a 
Kondo impurity and a symmetric Friedel impurity independent of any parameter. The 
only requirement is that the impurity has a resonance exactly at the Fermi level and 
that the level spacing is smaller than the resonance width. This result clarifies recent 
fidelity calculations. 

PACS: 75.20.Hr, 71.23.An, 71.27.+a , 05.30.-d 

In the late 1960s Anderson [1] showed that the potential of a weak impurity in a metal host 
changes the total n-electron wave function of the conduction electrons dramatically. Actually 
with increasing number N c of electron states (which is achieved by increasing the volume) 
the scalar product between the wave functions of the n-electron host without and with the 
impurity approaches zero. This phenomenon is generally called the Anderson orthogonality 
catastrophe (AOC). In recent years this phenomenon has been somewhat generalized and 
decorated with the romantic name fidelity. The generalization is that one applies the AOC 
to an arbitrary system which depends on one or several parameters A. If the system consists 
of electrons then it is described by its Hamiltonian. The Hamiltonian may contain a term 
which is proportional to a parameter A. Suppose that one can calculate the ground state of 
the system for A = and for finite A. Then the scalar product of the two wave functions 
is defined as the fidelity F^ c (0, A) of the system. Here N c is the number of conduction 
electrons states, which is proportional to the volume. The fidelity depends on the size of the 
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system and of particular interest is the limit for N c increasing towards infinity. If F^ c (0, A) 
approaches zero in this limit (the thermodynamic limit), then one faces an AOC. 

Our group studied recently the fidelity of the Friedel-Anderson impurity. This is an 
electron system with a d-atom as impurity. The energy of the d-electron lies at Ed below the 
Fermi level. If one removes a d-electron, i.e. creates a d-hole, then the conduction electrons 
can hop into the empty d-state with a hopping matrix element V s d- The d-hole possesses a 
finite life time before before it is refilled. Due to Heisenberg's uncertainty principle this 
life time broadens the energy level of the d-electron and transforms it into a d-resonance 
with a resonance width which is of the order of h/r^. In general the properties of a dissolved 
d-atom are more complicated because the different d-electrons repel each other due to the 
Coulomb interaction. In the theoretical investigation of such an impurity one studies (most 
of the time) a simplified model of a d-impurity with only two d-states, a spin-up and a 
spin-down d-state. Such an impurity was first studied by Friedel [2] and Anderson [3] and 
I call it the Friedel-Anderson (FA) impurity. The strength of the Coulomb interaction U 
represents a parameter A as introduced in the fidelity. For U = the impurity properties 
are much simpler. The impurity has a Friedel resonance at the d-energy Ed in each spin 
sub-band and is called a Friedel impurity. 

While the wave function of the FA impurity is quite complex the density of states is 
simpler and qualitatively sketched in Fig.l together with the density of states of a Friedel 
impurity. For both impurities the symmetric case is shown. In the Friedel impurity the 
resonance is positioned at the Fermi level (Ed = 0). In the FA impurity the d-state energy is 
positioned at Ed = —U/2 so that Ed and (Ed + U) lie symmetrically to the Fermi level. As 
a consequence there is an electron-hole symmetry. One obtains two d-resonances at roughly 
the energies Ed = —U/2 and Ed + U = +U/2 for spin up and down. These are known as 
the Hubbard resonances, and their width is twice the width of a Friedel resonance with the 
same s-d-hopping [I], [5], [6]. In addition one obtains a narrow resonance at the Fermi level 
which is generally called the Kondo resonance (see for example ref. [7]). 
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Fig.l: a) The electron band of a symmetric Friedel impurity with 
the d-resonance at the Fermi level. 

b) The effective band (density of states) after turning on the Coulomb 
repulsion U and setting Ed = —U /2. The Friedel resonances are 
transformed into broad Hubbard resonances, roughly positioned at 
Ed — — U/2 and E^ + U = +U/2, far away from the Fermi level 
and an extremely narrow Kondo resonance at the Fermi level. 



Two studies of the fidelity of the FA impurity have been published recently, one by 
Weichselbaum et al. [8] and one by our group [9]. In our investigation we calculated the 
fidelity between a symmetric Friedel impurity (Ed = 0, U = 0) and a symmetric FA impurity 
with finite Coulomb repulsion U and Ed = —U/2. When level spacing SE is smaller than 
the width of the Kondo resonance in Fig. lb then the fidelity did not change any more with 
increasing N. This raised a puzzled question by one referee of our paper: 'Does this mean 
that the simple Friedel resonance and the complex FA resonance have the same wave function 
(close to the Fermi level)'. In a way this question stimulated the present investigation. 

In the fidelity calculations one has on one hand to increase the number of electron states 
dramatically. On the other hand one needs to keep the number of states relatively small 
because otherwise the numerics requires an unacceptable computer time. These opposing 
requirements are optimally fulfilled by an ingenious trick applied by Wilson. One considers 
a system with 2 N I 2 electron states. The conduction band is half filled and symmetric to the 
Fermi level. For simplicity one assumes a constant density of states and divides all energies 
by the Fermi energy. Then the conduction band extends from (—1 : +1) as shown in Fig. 2. 
In the next step one divides the lower (and upper) half of the band geometrically into cells 
with decreasing width so that one obtains an energy frame. In Fig. 2 this energy frame has 
the values -1, -1/2, -1/4, -1/8, -1/16, -1/32, 0, 1/32, 1/16, 1/8, 1/4, 1/2, 1 and defines N 
energy cells € u (in Fig. 2 we have N = 12). The number N is always even because there are 



3 



as many negative as positive Wilson states. Each sub-division adds two states so that the 
number of sub-divisions or iterations is equal to N/2 . 

The four energy cells close to the Fermi level (in Fig. 2 (£ 5 , (£ 6 , (£ 7 , £ 8 ) may be considered 
to possess just one electron state each. The number of states per cell doubles with each 
step away from the Fermi level. (In Fig. 2 €\ and €12 have 16 electron states and the whole 
band has 64). Wilson reshuffled these electron states in each cell in such a way that one 
of them, the state cj, in the cell € u , carried the full interaction with the d-impurity. The 
remaining states in the cell have zero interaction and are neglected. The width of the four 
smallest energy cells in the vicinity of the Fermi level determines the effective energy spacing 
5E = 2-W 2 - 1 ) (in units of the Fermi energy). The effective number of electron states N e ff 
is determined by the energy spacing at the Fermi level and is given by N e ff = 2/6E=2 N / 2 
(for details see the discussion in ref. [9]). According to the construction the number N of 
Wilson states in FAIR is always even. 

The Wilson states have two great advantages: (i) One can represent a band of 2 N I 2 
electrons by just iV Wilson states and (ii) One can double the number of electrons by 
subdividing the cell directly below and the cell directly above the Fermi level into two equal 
halves. So the Wilson spectrum achieves the trick at doubling the effective number of states 
by adding just two states. 



ts 



+1 



+ 1/2 

+1/4 
+1/8 


-1/8 
-1/4 

-1/2 



-1 



Fig. 2: The Wilson description of a conduction band. 
The energies are normalized with the Fermi energy. 
The density of states is constant. The lower and upper 
half are subdivided into energy cells whose size is 
shrinks by A = 2 towards the Fermi level (at zero 
energy) . 



In this investigation I ask the question: Suppose that for a half-filled band with (N — 2) 
Wilson states we calculate the ground state of a given impurity and then repeat 

the calculation for a system with iV Wilson states for the same physical parameters, again 
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half filled. The occupation of each spin band is either (N — 2) /2 or N/2 which I define 
as n = N/2. The ground state with N Wilson states has one additional spin-up electron 
and one additional spin-down electron. The ground state can be easily expressed in 

the new basis. If we denote for a moment the two energy cells next to the Fermi level in 
the (N — 2)-basis as £_ and £ + then these cells are each split into half in the iV-basis. 
Consequently the amplitude in each of the new states is just l/y/2 of the amplitude in the 
original basis state. So the state is exactly transferred into the iV-state Wilson basis. 

The question in this paper is the following: What is the relation between *f> and \I/ jv i n 
the range of N where the fidelity is constant? I arrive at the following result: For sufficiently 
large N there is a single electron state such that ^ n — <?{<?{ ^ n-2 with remarkable accuracy. 
This conclusion will be derived in detail below. (In this paper I denote single-electron states 
such as \q) by their creation operator q^). 

The Kondo and the FA-impurity problem have been exactly solved with the Bethe ansatz 
|1(J] . [TTj . I am told that it is very hard to extract the wave function from the Bethe-ansatz. 
The presently most frequently used numerical method for the investigation of the FA impurity 
is the numerical renormalization group (NRG) theory, which was developed by Wilson [12] 
35 years ago and first applied to the FA impurity by Krishna-murthy et al. [13] in 1980. It 
derives the ground state through a large number of renormalization steps (of the order of 
50 to 100 steps). In each step the number of Slater states is increased by a factor of 16. (A 
Slater state is defined as the product of n single electron states). This yields a huge number 
of Slater states for the ground state and only a small number of the order of a few thousand 
Slater states are finally included in the calculation performed. 

For the actual calculation I use the FAIR ground state for the different impurities [H] . 
[T5] . [IB]. The FAIR technique has been developed during the past few years by the author. 
FAIR stands for Friedel artificially inserted resonance. The FAIR ground state represents 
a very good approximation for the Friedel-Anderson and the Kondo impurities. It has 
reproduced a number of numerical results with good accuracy and produced a number of 
new results, such as the polarization of the Kondo cloud [T7] , oscillations in the Kondo cloud 
|18j . Friedel oscillations of the FA impurity [TS], and, in the magnetic pseudo-ground state, 
the magnetic moment [2D] which roughly corresponds to the mean field result with half the 
Coulomb energy U because of the reduced density of states in the d- resonances [6j. A review 
of the FAIR theory is given in ref. [21]. 

The singlet ground state of the FA impurity, which consists of eight Slater states, was 
used to calculate the fidelity between the symmetric Friedel and the symmetric FA impurity 
[9]. Its wave function is essentially the superposition of two magnetic states ^ms with 
opposite magnetic moment. The magnetic state (with net spin up) has the form 

OfCW (1) 

where aL and feL are two artificial resonance states or fair states in the spin-up and down 
bands. It is defined and investigated in ref. [20] and [6J within FAIR formalism. Its magnetic 
moment is \B 2 — C 2 \ in units of (j,b- 
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The composition, for example, of a ^ in terms of the N Wilson states c u ^ is 

TV 
v=l 

The fair states assume the effective interaction with the impurity. Since the electron system 
has two spin sub-bands one needs two fair states, the state for the spin-up sub-band and 

the state 6^ for the spin-down sub-band. The remaining states aL in the spin-up sub-band 

are made ortho- normal to each other and to a^. In addition their free electron Hamiltonian 

matrix (a^Q \ H° \ a j^fi ^ is sub-diagonalized (excluding the row and column with matrix 

elements). As a result the state aL becomes an artificial Friedel resonance. The fair state 

aL (and bL) determines uniquely the remaining band states aL (and oj,) for % > 1 which 
form a new conduction band with one less electron (for each spin). The half-occupied bands 
are represented by 

n— 1 n— 1 

lOato^Hn^n^ 

i=l i=l 

where n = N/2. 

The FA ground state is the normalized superposition of the magnetic states with net spin 
up and net spin down. 

^SS = *AfS t + ^msi (2) 
The state ^msi is obtained by reversing all spins in equ. ([I]) (the spins are ordered in the 
same fashion as in equ. (JT])). The main numerical task is to find the optimal fair states ag 
and £»q (which occur now in both spin directions). When this is achieved by variation the 
full wave function can be easily constructed. 

The FAIR technique has a number of advantages: (i) Two single electron states, the 
fair states and b\ determine the full bases of the electron bands parallel and anti-parallel 
to the impurity spin. Each fair state requires only a small number (of the order of 40) of 
coefficients of Wilson states, (ii) The ground states for the FA and the Kondo impurity 
consist of a small number of Slater states, (iii) The d-state occupations in the different 
Slater states is well separated insofar as each Slater state possesses either zero, one, or two 
d-electrons. 

So far the error margins of the FAIR technique have not been quantified. However, the 
quality of reproducing previous results justifies the use of this rather transparent method to 
predict new phenomena and uncover relationships and coherences which were not transparent 
before. 

In this paper I denote the singlet ground state \& ss m the basis of iV Wilson states as 
^n- For the FA-impurity ^ n can be written as a sum of (eight) Slater states with N/2 
spin- up and N/2 spin-down states. Each Slater state S N is the product of N/2 spin- up and 
N/2 spin-down electron states (creation operators applied to the vacuum state fi) and has 
a coefficient as N - I denote the sum of all these Slater states including their coefficients as 
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A^Q = ^jv- Similarly one can express the ground state for the (N — 2)- Wilson basis as 
V N _ 2 = B j N _ 2 Q. Then one can multiply the state with the adjoint operator of B^ N _ 2 ,i.e 
the corresponding annihilation operators and form the state 



B ] 

- D iV-2 



f A ] N n = B N _ 2 d N tt = a Q Qln (3) 



This procedure yields a two-electron state chqQIQ with one electron in each spin, (ocq is the 
amplitude and Q\ is normalized). In the final step I try to express the two-electron state 
Ql as the product of a spin-up and a spin-down single-electron state ql and q\ with identical 
orbits, 

a Q Q\Q = a 2 q q\qp (4) 

The value ot q is the amplitude of the (normalized) state q\. The absolute value \a q \ is less 
than or equal to one. This can be seen in the following way: We form the iV-particle state 
B^ n _ 2 Q\Q and take the scalar product between B^ n _ 2 Q\Q and A^ N Q = ^at- 

Bl_ 2 Qitt\A]<ri) = (Qln\B N - 2 A* N n) = (Qln\a Q Qln) = a Q 



where equ.© is used. The scalar product is equal to «q. The absolute value of the scalar 
product between two normalized iV-electron states is less than or equal to one. Therefore 
one has \olq\ < 1. Under optimal conditions the two-electron state Q\ can be factored into 

two single electron states, i.e. aqQ\ — [oiqQ^ (^qlfj ^- O ne expects that the coefficient 
\a q \ is also less than or equal to one. 

After the two ground states \1/ jy and \I/ at__ 2 are constructed in form of eight Slater states 
the calculation of B^_ 2 A^ N VL consists mainly of scalar products between different 



electron states (n = N/2). The latter are determinants of single-electron scalar products. 
The resulting coefficients of oqQ^ form a quadratic (N + 1) x (JV + l)-matrix in terms of 
A" Wilson states plus one d-state for each spin. In the last step the two-electron state oiqQ\VL 
is split into the product (a q q^j (otqq^J This procedure is remarkably easy. Already the 
square root of the diagonal elements in o.qQ\VL yields the absolute value of the coefficients 
of (ct q qf^ ^, and the sign follows from the non-diagonal elements. 

Table I shows the value \a q \ for a number of systems. The first column gives the impurity, 
the second the s-d-hopping strength, and the third column the number A" of Wilson states. 
There are three different values of A" used, A" = 42, 32 and 22. The corresponding values of 
the smallest energy spacing are 2~ 20 ps 1CT 6 , 2~ 15 ps 3. x 10~ 5 and 2 -10 ps 1CT 3 . 

For the FA impurities I used U — 1 and = —0.5. In all these cases the split of the 
two-particle state ocqQ^Q into two identical single particle states with opposite spin worked 
very well. The coefficient a q in column four is very close to one. The fifth column gives 

the Gaussian deviation E rr between aqQ 2 Q, and (a q q^\ (a q q^J fi. Obviously as soon as the 

smallest energy in the Wilson spectrum SE is sufficiently small compared to the resonance 
width of the impurity at the Fermi level then the state q 1 is well defined and almost perfectly 
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normalized. This means that the state can be constructed from the state by the 

relation 

V N q\q\^ N - 2 (5) 



impurity 


|v sd | 2 


N 


\a q \ 


E„ 




c g (0) 


Friedel 


10" 4 


42 


0.995 


io- 13 


0.814 


-0.528 


Friedel 


0.03 


22 


0.995 


io- i3 


0.813 


-0.529 


FA 


0.03 


42 


0.990 


3 x 10 4 


0.814 


-0.525 


FA 


0.04 


42 


0.998 


5 x 10~ 5 


0.814 


-0.528 


FA 


0.05 


22 


0.979 


2 x 10 3 


0.812 


-0.504 


FA 


0.05 


32 


0.997 


2 x 10~ 4 


0.814 


-0.527 


FA 


0.05 


42 


0.998 


2 x 10~ 5 


0.812 


-0.533 


Kondo 


J = 0.10 


42 


0.992 


2 x 10~ 4 


0.8062 


-0.533 



Table I: The first column gives the impurity, the second the 
s-d-hopping strength, the third column the number N of Wilson 
states. The fourth column gives the coefficient \a q \ of the 
extracted state . The next column gives the Gaussian 
deviation _E rr .The following two columns are explained in the text. 



In all cases the extracted state q* is mainly composed of basis states close the Fermi level. 
In table I the two last columns give the amplitudes of the state q* below and at the new 
Fermi level in terms of Wilson states. These amplitudes agree almost perfectly. 

Fig. 3 shows the coefficients of q 1 with respect to the Fermi level for pure Friedel impurities 
with N = 42 and very different s-d-hopping strengths \V s d\ ■ It also shows the corresponding 
coefficients of q* for FA impurities with N = 42 and different \V s d\ 2 (which yield many orders 
of magnitude different Kondo energies). In addition Fig. 3 shows the effect of different iV 
for a FA impurity with \V sd \ = 0.05, U = 1 and E d = -0.5. The Kondo energy is about 
3.9 x 10~ 3 . For iV = 22 the smallest energy spacing is about 10 3 . Still the N = 22 curve 
is relatively close to the universal curve although table I shows that its \a q \ = 0.979 is less 
close to one than for larger N. For values of iV where the smallest level spacing is of the 
order of the Kondo energy or larger the value of a q decreases. Then the state ^ N can no 
longer be expressed as q\q\^ n-2 and the splitting of Q\ into two single electron states with 
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opposite spin becomes meaningless. 




Fig.3: This figure shows the coefficients of the state for 
three different impurities, F =Friedel, FA=Friedel-Anderson 
and K=Kondo. The coefficients are counted from the Fermi 
level. Most curves are for N = 42, but one curve is for a 
FA impurity with \V ad \ = 0.05 and N = 22. The curve shown 
is super-universal because it is independent of the kind of 
impurity (F,FA,K), independent of N, and independent of \V s d 
and other parameters as long as there is a resonance exactly 
at the Fermi level. 



I also performed the same calculation for a Kondo impurity. The parameters are collected 
in table I as well. The same behavior is observed. Fig.3 includes the coefficients of a Kondo 
impurity with J = 0.10 and N = 42 as open stars. 

The interpretation of the above results is that the complexity of the solutions of the 
FA and the Kondo impurity is not felt at energies sufficiently below the Kondo energy. 
We know from Wilson's renormalization calculation that the structure of the ground state 
changes dramatically when the smallest level spacing 5E is of the order of the Kondo energy. 
After this metamorphosis the renormalization approaches a fixed point. From this behavior 
Nozieres (22] developed the Fermi liquid theory of the Kondo impurity (which applies also 
to the FA impurity). The present calculation provides quantitatively a universal state q\ 
which has to be incorporated into the ground state for both spins when the Wilson basis is 
increased by 2 states. The fact that this is the same state for any impurity with a resonance 
at the Fermi level demonstrates that the Kondo and FA impurities have a resonance at 
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the Fermi level. The result suggests an induction rule for the ground state with increasing 
number of Wilson states. 

The accuracy of relation (J5J) is given by the scalar product between the left and right side. 
This is for the FA impurity and N = 42 already 0.996 (i.e equal to \<y q \ 2 ) and approaches 
the value one with increasing N. 

The fact that is the same state for any impurity with a resonance at the Fermi level is 
at the heart of our earlier discussed fidelity calculations. The fidelity is the scalar product 
between the ground states of the symmetric Friedel and FA impurities, which 
is compared with the corresponding scalar product in the (N — 2)-basis \&n-2\^n-2/- (J n 
both cases the bands are half-filled). According to our present result we can express the 
latter as 

This relation confirms that the fidelity between the symmetric FA and the symmetric Friedel 
impurity approaches a finite value with increasing N and does not experience an Anderson 
orthogonality catastrophe. The equation (j3J) yields an asymptotic recursion formula to con- 
struct from a given ty^o where No is sufficiently large so that the smallest level spacing 
in the Nq- Wilson basis is smaller than the Kondo energy. 



Conclusion: In the Wilson basis one can easily introduce two new states by splitting each 
of the states directly below and above the center of the band into two states of equal weight, 
creating out of two old states four new ones. This reduces the smallest energy spacing 
by a factor of two and increases the effective number of states N e ff by a factor of two. 
The relation between between the ground states and \I/jv of the two half filled bands 

is investigated for several symmetric impurities, the Friedel, the Friedel-Anderson and the 
Kondo impurity. The ground states are calculated with the FAIR technique. We observe 
that for sufficiently large N the ground state can be obtained from by multiplying 

by two single particle states ql and q\. The composition of these states g^with respect 
to the Fermi level is independent of N (for sufficient large N) and is independent of the 
resonance width of the impurity. This qualifies the state q^ to be called universal. The state 
q^ is not just universal for a given kind of impurity. It is super-universal because it is even 
independent of the kind of impurity as long as the impurity has a resonance directly at the 
Fermi level. In particular it is independent of the complexity of the ground state of the FA 
or the Kondo impurity. All the complexity happens further away from the Fermi level. This 
result gives a simple confirmation and explanation of our previous fidelity calculations [9]. 
A similar investigation for the asymmetric FA impurity is in progress. 
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